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Abstract 

The relativistic field theory model of the deuteron suggested by us previously is 
applied to the calculation of the reaction rate of the low-energy two-proton fusion 
p + p^D + e"*" + fe. The theoretical prediction of the reaction rate obtained 
is 2.9 times larger than given by the potential approach. This leads to a strong 
suppression of the high energy solar neutrino fluxes. 
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1 Introduction 



The solar neutrino problem has its origin in the strong discrepancy between predictions 
of the Standard Solar Model (SSM) [1,2] and the observed solar neutrino fluxes [2]. The 
experimental solar neutrino fluxes turned out substantially less than that predicted by 

SSM [1-3]. 

One of the solution of this problem is to lower the temperature in the center of the 
Sun in comparison to that predicted by the SSM: Tc = 1.55 x 10^ K [1,2]. Indeed, due 
to strong dependence of the solar neutrino fluxes on Tc just a 20-^30% diminishing of 
leads to a decrease of the neutrino fluxes by more than in order of magnitude [3]. 

In order to reduce one can resort to the change of physical and chemical phenomeno- 
logical inputs which determine the structure of the star [3]. We will not enumerate here all 
of them, and relegate the reader to Ref . [3] . We just pick up the temperature associated 
with the magnitude of the reaction rate of the low-energy two-proton fusion, i.e., p + p 
^ D + c+ + z/c. Since the solar luminosity should be kept constant, an enhancement of 
the reaction rate magnitude of the two-proton fusion implies a lower T^. 

The low-energy electroweak model and non-relativistic approaches, having been ap- 
plied to the computation of the contribution of strong interactions to the matrix element 
of the p + p — > D + W"*" transition does not leave a room for the substantial change of 
the cross section magnitude [1,2,4,5]. As has been noted in Ref. [3] they differ from the 
mean value by no more than 3 %. 

Recently we have suggested a phenomenological relativistic held theory model of the 
deuteron [6,7] . In the relativistic field theory model of the deuteron the strong interactions 
of particles are described by one-nucleon loop diagrams. In Ref. [7] we have computed 
the cross section of the radiative proton-neutron capture n + p — 5> D + 7 for thermal 
neutrons and the cross section of the two-proton fusion p + p ^ D + e+ + z/g. In the case 
of thermal neutron radiative capture our result agrees well with the experimental data 
and the non-relativistic potential approaches. However, in the case of the two-proton 
fusion we have found an enhancement of the reaction rate by a factor of 2.9 with respect 
to the non-relativistic potential approaches [1,2]. 

This result is due to our model approach using one-nucleon loop diagrams for the 
description of strong low-energy interactions of the deuteron to other particles. It is 
well-known that such fermion-loop diagrams should possess anomalies [8-12]. As usually, 
these anomalies dominate for the amplitudes of strong low-energy interactions of hadrons 
[8,13-15]. In the case of the two-proton fusion we encounter the dominance of the anomaly 
of the AAV one-nucleon loop diagrams, that is the diagrams having two axial-vector and 
one vector vertices. Such an anomalous contribution produced by vacuum fluctuations 
of virtual nucleons has completely a quantum field theory nature and cannot be restored 
within the potential approach describing strong low-energy interactions in terms of the 
overlap integral of the wave-functions of the protons and the deuteron. The result of the 
computation of the AAV-anomaly depends on the shift of the virtual nucleon momentum. 
This is similar to the well-known Adler-Bell-Jackiw-Bardeen anomaly [8,9] described by 
the AVV-fermion diagram. This ambiguity can be fixed by the requirement of gauge 
invariance under gauge transformations of the deuteron field. This is similar to the removal 
of the ambiguity appearing for the computation of the AVV-anomaly. 

The paper is organized as follows. In Sect. 1 we adduce the computation of the cross 
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section and the reaction rate of the two-proton fusion p + p — D + e+ + z/q [7]. In Sect. 2 
we apply the obtained cross section to the discussion of the solar neutrino problem. In 
the conclusion we discuss the obtained results. 



2 Cross section of two— proton fusion 

In the relativistic field theory model of the deuteron the effective Lagrangian describing 
strong interactions of the deuteron with nucleons reads [6,7] 

Aot(x) = -\dI^{x)D^^{x) + MlDl{x)D^{x) - igy[p{x)rn%x) - 

-n{x)j^p''{x)]Dfj_{x) - igy[p''{x)'y^n{x) - n'=(x)7^p(x)]Dj,(x) + (2.1) 
+p{x){i'^^d^ - Mn)p(x) + n{x){i'^^d^ - M^)n{x), 

where D^{x), p{x) and n{x) are the interpolating fields of the deuteron, proton and 
neutron, respectively, D^^{x) = d^D^{x) — d^D^{x), = 2Mn — is the mass of the 
deuteron, Mn = 938 MeV is the mass of the proton and neutron, and Ed is the binding 
energy. In our approach is given by [6,7] 

The phenomenological coupling constant gy is connected with the electric quadrupole 
moment of the deuteron Qd by the relation 

By using the experimental value (QD)exp = 0.286 fm^ [16] and Mq — 1876 MeV one can 
estimate the value of gy with gy = ±11.3. Without loss of generality one can use the 
positive sign, i.e., gy = 11.3 [7]. The quantity Ad = 64.843 MeV is the 3-dimensional 
cut-off in one-nucleon loop diagrams. The magnitude of Ad is obtained by the fit of the 
binding energy e-^ = 2.225 MeV[16]. We identify 1/Ad with the effective radius of the 
deuteron [6,7]. At Ad = 64.843 MeV we get td = 1/Ad = 3.043 fm [7], which agrees 
well with the average value of the deuteron radius, i.e., < r > = 3.140 fm [17]. 

First, for the computation of the effective Lagrangian describing the low-energy pp- 
scattering we perform keeping the one-pion-exchange contribution 

>C^^(x) = - iPixh'pi^)] ipi^h'pix)] . (2.4) 



Here we have denoted g-n-NNiv) = gTTNNC{v) and C{v) = ^^j2'Ka/vexp{—na/v), f is a 
relative velocity of interacting protons, which takes into account the Coulomb repulsion 
between protons at low energies [4]. 

Applying the Fierz transformation we bring the Lagrangian (|2.4|) to the form 
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>^e&(^) = %^ {\p{x)f{x)mx)p{x)] + \p{xh'f{x)mxh'p{x)] + 

1 , ^2.5) 

+ b(^)7M7V(a;)]b'(a^)7Vp(a;)] + ^\p{x)a^,^p''{x)]\p''{x)a'"'p{x)] }. 

The process p + p— ^D + e^ + i^e should run via the intermediate W-boson exchange, 
i.e., p + p ^ D + — > D + e"*" + z^e- The Lagrangian describing the electroweak 
interactions of the W-boson with proton, neutron, positron and neutrino reads [18,19] 

C(^) = - 9A7'Hx) + u,{x)Y{l - l')e{x)] W,{x) + h.c. (2.6) 

Here is the electroweak coupling constant connected with the Fermi constant Gp — 
1.166-5 GeV"^ and the W-boson mass Mw by the relation [18,19] 



9^ G 



F 



8M2, V2' 



(2.7) 



where gp^ = 1.260 ± 0.012 is the axial-vector coupling constant [16] describing the renor- 
malization of the weak axial-vector hadron current by strong interactions. 

The interaction \p{x)'j^'y^p'^{x)]\p'^{x)'y^'y^p{x)] gives in our approach the main contri- 
bution to the amplitude of the transition p + p — > D + W. The corresponding one-nucleon 
loop diagrams are depicted in Fig. 1. 

The effective Lagrangian defined by the diagrams in Fig. 1 is given by 



/ 



d^x Cpig,i{x) = 

I / ^-J0^ m^halM^)] ^i(^l) W!{X.) X (2.8) 

-iki-xi -ik2-X2 ^i{ki+k2)-x ^ „ 9W 5'^Nn('^) 9V 



where 



J''''\k^M\Q) = (2.9) 



1 1 1 



M^-k-Q M^-k-Q-ki - k - Q - ki - k2 . 



The momentum integral J°''^''{ki, k2] Q) defines the structure of the effective Lagrangian, 
the 4-vector Q — aki + bk2 describes an arbitrary shift of the virtual momentum of 
the nucleons in the diagram in Fig. 1. The parameters a and b are free. It should be 
emphasized that a and b are not phenomenological parameters of the model, they appear 
only due to the properties of fermion-loop diagrams under the shift of virtual momenta 
[8,11,12]. 

Keeping only the anomalous part of the momentum integral [13-15] and the leading 
order of the expansion in the powers of the momentum k2 coinciding with the momentum 
of the leptonic pair that is soft, removing then ambiguities connected with the shift of the 
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virtual momentum by the requirement of the gauge invariance under gauge transforma- 
tions of the deuteron field, we arrive at the structure function [7] 



J^^^{k,, h- Q) = 3 {k'^g''^ - k^.gn- (2-10) 
The effective Lagrangian defined by the structure function ( |2.10| ) is given by [7] 

/:F.g.i(x) = - 1^ |j ^i^^^^ mx)ri'p{x)]. (2.11) 

In turn the effective Lagrangian describing the low-energy process p + p^D + e'^ + Ug 
reads 

^efr(x) =9A^^ W^^'"^"^^ ^i'^^'^^ mx)Yl'p{x)], (2.12) 

where j^{x) = i>e(a;)7^(l — 7^)e(x) is the leptonic electroweak current. 

In the low-energy limit when the 3-momenta of protons tends to zero the effective 
Lagrangian (|2.12|) can be reduced to the expression 

£efr(x) = ^gAM^^ ^ ^||^ f{x) DI{x) mx)j'p{x)] , (2.13) 
where we have used the relations 

[p^{x)Yi'p{x)] ^ -g'^lp'ixh'pix)] 

Dl,{x)^-iM^Dl{x) ^ ■ > 

that are valid in the limit of low 3-momenta of protons and the deuteron. 
The amplitude defined by the effective Lagrangian ( |2.13D reads 



A^(p + p ^ D + e+z/e) = iC{v) gAM^ 



Gf ^gv 



NN 



y/2 47r2 8M2 

X eliQ) [u{k,)r{l - mPi)lMP2)l (2.15) 

where u{ky) and f (A^e) are the Dirac bispinors of the neutrino and positron. For conve- 
nience we have separated the factor C{v) describing the Coulomb repulsion. 

The amplitude Eq. ( |2.15| ) takes into account only the contribution of the one-pion 
exchange. Since at low energies the pp-system is in the ^Sq state with the isospin equal 
T = 1, there is the contribution of the pole on the unphysical sheet. Following the com- 
putation of the cross section of the radiative neutron-proton capture [7] and normalizing 
the contribution of the pole on the unphysical sheet to the ^Sq pp-scattering length as 
we get 

A^(p + p ^ D + e+z/e) = I C{v) gAM^% ^ e* (Q) [u{k,)i^{l - i'')v{k,)] x 



47r 



as 
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The experimental value of as is given by as = (—17.1 ± 0.2) fm [16]. 

For the correct description of the strong low-energy pp coupling to the deuteron 
the amplitude Eq. (|2.16|) should be unitarized [7]. In our approach we can perform the 
unitarization by summing up an infinite series of proton-proton bubble diagrams [7] . This 
gives 

^(p + p ^ D + e+z/e) = t C{v) QaM^^ ^ e*(Q) [uiKh^il - ^'Hh)] x 



mi 



8V2nMl as 



X 



8^27^ Ml as \ Ml 



^/2 47r2 

— [m^(Pi)7^w(P2)], 



(2.17) 



where v is the relative velocity of the protons [7]. 

The cross section of the low-energy p + p^D + e+ + z/e scattering is defined 



(p + p^D + e+ + z/e) = -TT^ / |A^(p + p^D + e+ + z/e)|2 X 



(27r)32ED (27r)32Ee {2iif2E^ 



(2.18) 



where |A^(p + p — D + e+ + z/e)p is the squared amplitude averaged over polarizations 
of protons and summed over polarizations of final particles [7] 



|A<(p + P - D + e+z.,)|2 = C\v) glM^^^^ 

ion 



X 





5^NN 
_8M2 


- 2 


/ 8V27r Ml as 


r 


1 + 


n2 1 

i/TrNN 

8M2_ 


■( 


^ 8V2n Ml as ^ 


47r2 



QaQl3 



X 



X tr{(^e - me)7"(l - 7')^.7^(1 - 7')} ^ H{p, - Mn)7'(P2 + Mn)7'} 



c . 3 e . 



X 





[^2 1 
i/vrNN 

_8M2_ 


8V2nMl as y 


1 + 


[^2 1 
i/TrNN 

_8M2_ 


2 


( 8V27r Ml as V M^ , 
[ gl^^ M^J 4vr2^ 



(2.19) 



Now we should carry out the integration over the phase volume of the final D e^i^e^state 



d'Q 



d kn 



d k]. 



(27r)32ED (27r)32^e {27t)^2E^ 
x{27r)U^'\p,+P2 -Q-K-K) (e,E, - ■ k^j 

\— r ^ El - mlE,{e^ - dE, 



327r2 



9607^2M^ 



/(O, (2.20) 
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where ^ = m^jsYy and the function /(^) is defined by the integral 



30 



/' 4^^^x (1 - = (1 - ^ - 4a 0^ + 



= 0.776 



(2.21) 



and normahzed to unity at ^ = 0, i.e., /(O) = 1. 

The cross section of the low-energy two-proton fusion p + p — > D + e"*" + i^e is given 

by 



(t(p + p ^ D + e+ + i/e) = 



X 



8MI 



V 12807r5 

' ' 8V2nM^ as V 



f 



X 



1 - 



Ml 



N , 



1 + 



SMI 
1 



'^_8V27rMlas_yM^^2 



1 + 52743 



X 



X 1.34 X 10~^^cm^ 



(2.22) 



The cross section is calculated in units of h = c = 1. The appearance of the factor C^(i') 
agrees with the result obtained by Bcthc and Critchfield [4]. 

The reaction rate of the process p + p— >D + e"'" + t'e reads [7] 



< V ct(p + p ^ D + e+ + z/e) >= 1.03 x 10 



-39 



1 2 
X — X — X 

a TT 



X 



X 



1 + 52743 



2^2 



1.03 X 10-=^^ gg cm^s-\ 



where r is connected with the temperature T [4] 



N 



1/3 



4kT 



1 + — 

T 



X 



(2.23) 



(2.24) 



and k = 8.62 x 10^^^ MeV ■ is the Boltzmann constant. Following [4] we have com- 
puted the integral over relative velocities of the protons by a saddle point approximation. 

Setting T — Tc — 15.5, measured in unites of 10^ K, where Tc is the temperature of 
the solar core in the Standard Solar model [1,2], we obtain r = 13.56 and the reaction 
rate is given by 



< v (t(p + p ^ D + e+ + i/e) > = 3.44 x 10~^^ cm^ s 



3 o-l 



(2.25) 



This value differs by a factor of 2.90±0.87 from that calculated within the potential 
approach [1,2] (see also Ref. [5]). In order to reconcile our result with the solar luminosity 
we must assume that the temperature in the solar core equals Tc — 13.8+oi ^ 10^ K [3] 
here we assumed. As has been remarked in Ref. [3] the increase of the magnitude of the 



7 



cross section of the two-proton fusion leads to a decrease of the temperature in the solar 
core, and provides a strong suppression of the solar neutrino flux. 

The enhancement of the amplitude of the p + p — > D + W transition found in our 
approach is related to the computation of the amplitude in terms of one-nucleon loop 
diagrams. Indeed, the structure function Eq. (|2.10|) defining the effective Lagrangian of 
the p + p — s> D + W transition is due to the contribution of the anomalous part of the 
AAV one-nucleon loop diagram [11,13,15]. Such an anomalous contribution, produced by 
vacuum fluctuations of virtual nucleons, has a quantum fleld theory nature and cannot 
be obtained within the potential approach describing strong low-energy interactions of 
the protons and the deuteron in terms of the overlap integral of the wave-functions of 
two-protons and the deuteron. The ambiguity of the computation of the AAV-anomaly, 
produced by the shift of virtual nucleon momentum, has been flxed by the requirement 
of gauge invariance under gauge transformations of the deuteron fleld [7]. This is very 
similar to the removal of the ambiguity appearing for the computation of the Adler-Bell- 
Jackiw-Bardeen anomaly [20]. 



3 Solar neutrino fluxes 

The solar neutrino signals Sci and Sca obtained in the Gallium [21] and Chlorine [22] 
experiments, respectively, are given by 

SGa(^^ > 0.233 MeV) = (77 ± 8) SNU, 

Sci(^z. > 0.814 MeV) = (2.55 ± 0.25) SNU, (3.1) 
SKam(^^ > 7 MeV) = (2.73 ± 0.38) ■ lO'' cm-^s-\ 

where the Solar Neutrino Unit (SNU) is deflned as one reaction per second per 10^^ 
atoms [3] . We have adduced also the data given by KAMIOKANDE experiment [23] (see 
also Ref. [3]). Due to the energy threshold in KAMIOKANDE experiment this neutrino 
detector is only sensitive to neutrinos. 

The solar neutrino fluxes where i = pp, pep, ''Be, ^^N, ^^O and ®B, etc. can be 
represented in the form of a power-law behavior (see Eqs. (96) in Ref. [3]), i.e., 

$. = (3.2) 
is the neutrino flux, predicted by SSM, and the parameter x in our deflnition reads 

X = <--iP + P^^ + -^ + -eh> ^ 2.90 ± 0.87, (3.3) 

< a(p + p -> D + e+ + i/e)s >* 

where < v cr(p + p D + e+ + z/c)s >* is the quantity calculated in the potential approach. 

The values of the parameters ai are given in Table 1 and can be found in Table X of 
Ref. [3]. In Table 1 we also adduce the neutrino fluxes that should give the contributions 
to the signals detected in Gallium and Chlorine experiments. We have normalized our 
predictions to the results obtained within SSM [2]. 

It is seen that our model explains reasonably well the experimental data of the Gallium 
experiment. For the neutrino flux measured in the Chlorine experiment our model leaves 
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room for possible other contributions. One can see that our prediction for SKam(-E'i^ > 
7 MeV) = ^b[W^ cm-2s-^] = 0.37;g;^? (see Table IX of Ref. [3]) is too small when 
compared with the KAMIOKANDE experimental data. 

4 Conclusion 

We resume the obtained results. We have shown that the relativistic field theory 
model of the deuteron, describing the strong interactions of the deuteron with hadrons 
in terms of one-nucleon loop diagrams, has provided contributions of strong interactions 
to the matrix element of two-proton fusion p + p— >D + e+ + z/e, that has led to 
an enhancement of the reaction rate by a factor 2.9 with respect to the computations 
performed within non-relativistic approaches. 

This explains reasonably well the experimental data of the Gallium experiment. For 
the neutrino flux measured in the Chlorine and KAMIOKANDE experiment we obtain 
smaller values. This implies that the enhancement of the cross section of the two-proton 
fusion obtained in our model leads to a suppression of the high energy neutrino fluxes. 

In the conclusion let us discuss the estimate of the theoretical uncertainty of the 
relativistic field theory model of the deuteron. Unfortunately, in our model there is 
not any small parameter like a = e^/An = 1/137, the fine structure constant in QED 
or a large parameter N, the number of quark colors in a multicolor extension of QCD 
with SU{3)c — > SU{N)c, that allows to apply the expansion in powers of so-called 
large N-expansion, to quark-gluon diagrams describing strong low-energy interactions of 
hadrons. 

Our effective approach is based on the one-nucleon loop diagram approximation for 
the description of both the self-interactions of the deuteron and coupling of the deuteron 
to other particles. By describing the self-interactions of the deuteron at the one-nucleon 
loop approximation we have fitted all parameters characterizing the physical deuteron. 
In this case the only way to estimate the theoretical uncertainty of the approach is to 
compared the one-nucleon loop data with the two-nucleon loop corrections. In Refs. [6,7] 
we have calculated the two-nucleon loop correction to the binding energy of the deuteron. 
This correction reads [6,7] 

The numerical value of the two-nucleon loop correction makes up about 16% of the binding 
energy e-o = 2.273 MeV calculated in the onc-nuclcon loop approximation [7]. Thus the 
magnitude A = ±16% might be accepted as a theoretical uncertainty of the relativistic 
field theory model of the deuteron for amplitudes. The theoretical uncertainty of the cross 
section should be about 30 %. This agrees well with our prediction for the cross section of 
the radiative neutron-proton capture. Indeed, the theoretical value (T(n + p— >D + 7) = 
276 mb [7] is in agreement with the experimental data a{n + p — * D + 7)exp = (334.2 ± 
0.5) mb with an accuracy of 22 %. By taking into account the theoretical uncertainty, the 
cross section of the radiative neutron-proton capture should read: a{n + p — >^ D + 7) — 
(276 ±83) mb. Taking into account the theoretical uncertainty of the model we predict for 
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the two-proton fusion reaction rate an enhancement by a factor of 2.90 ± 0.87 compared 
to the potential approach. 
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Table 1. Contributions from the main components of the neutrino fluxes to the signals 
(SNU) in the Gallium [21] and Chlorine [22] experiments according to the SSM [2] and our 
approach (see Table VII of Ref. [3]). The errors are due to the assumed 30% uncertainty 
of the reaction rate for p + p — > D + e"^ + i^e- The power-law parameters ai have been 
taken from Table X of Ref. [3]. 
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Figure Captions 

• Fig. 1 The contribution of \p{x)^^^^p'^{x)]\p^{x)^^^^p{x)] to the amphtude of the p 
+ p — > D + e+ + i^e scattering. 
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